Non-Abelian Geometric Quantum Memory with Atomic Ensemble 



Y. Li \ P. Zhang \ P. Zanardi 2 < 3 , and C. P. Sun l < a < b 
1 Institute of Theoretical Physics, The Chinese Academy of Science, Beijing, 100080, China 
2 Department of Mechanical Engineering, Massachusetts Institute of Technology, Cambridge, MA 02139 
3 Institute for Scientific Interchange (ISI), Viale Settimio Severo 65, 10133 Torino, Italy 

We study a quantum information storage scheme based on an atomic ensemble with near (also ex- 
act) three-photon resonance electromagnetically induced transparency (EIT). Each 4-level-atom is 
coupled to two classical control fields and a quantum probe field. Quantum information is adiabati- 
cally stored in the associated dark polariton manifold. An intrinsic non-trivial topological structure 
is discovered in our quantum memory implemented through the symmetric collective atomic excita- 
tions with a hidden SU(3) dynamical symmetry. By adiabatically changing the Rabi frequencies of 
two classical control fields, the quantum state can be retrieved up to a non-abelian holonomy and 
thus decoded from the final state in a purely geometric way. 
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Quantum information storage is a physical process to 
encode the state of a quantum system into the state of 
another system referred to as a quantum memory [1]. 
Compared to the original quantum system the quantum 
memory should possess a large decoherence time for effec- 
tive storing of quantum information. Moreover, the orig- 
inal state of the quantum system should be retrievable 
from the encoding quantum memory state. By means of 
quantum memory one can transport quantum informa- 
tion from place to place within the decoherence time. Re- 
cently ensemble of A-type atoms has been proposed [2-4] 
as a candidate for practical quantum memory. The idea 
is to store and transfer quantum information contained in 
photonic states by the collective atomic excitations. This 
approach is based on the phenomenon of electromagneti- 
cally induced transparency (EIT) [5] . Some experiments 
[6,7] have already demonstrated the central principle of 
this technique, namely, the reduction of the group veloc- 
ity of light. 

Most recently a system with quasi-spin wave collec- 
tive excitations of many A-type atoms fixed in "atomic 
crystal" has been considered as a candidate for a robust 
quantum memory [8]. A hidden dynamical symmetry of 
such a system is discovered and it is observed that in cer- 
tain cases [9] the quantum state can be retrieved up to 
a non-abelian Berry phase, i.e., a non-abelian holonomy 
[10-15]. This observation extends the concept of quan- 
tum information storage. Quantum information storage 
of photonic states with this topological character can be 
implemented in an atomic ensemble with off-resonance 
EIT. In such a case the stored state can be decoded 
in a purely geometric way. However, this non-abelian 
holonomy is in some sense trivial due to the fact that the 
quantum storage space splits into an orthogonal sum of 
invariant one dimensional subspaces. 

In this letter, we shall describe a quantum information 
storage protocol based on a truly non-abelian holonomy. 
To this aim we will consider an ensemble of N 4-level- 
atoms [13,14], where two meta-stable states are coupled 
to the excited state by two classical control fields respec- 



tively while the ground state is coupled to the excited 
state by a quantum probe field. In the large N limit with 
low excitation, a three-exciton system is formed by the 
symmetric collective excitations from the ground states 
up to the excited state plus the two virtual excitations 
from the two meta-stable states to the excited state. It 
is easy to prove that these three collective excitations in- 
deed behave as three bosons in the large N limit with 
low excitation. Intertwining between the excited state 
and two meta-stable ones, the collective operators gen- 
erate an SU(3) algebra. Based on the spectrum gener- 
ating algebra theory [16] associated with this SU(3), we 
construct the degenerate eigen-states of the three-mode 
exciton-photon system. In particular the collective mani- 
fold of dark states can be shown to split into dynamically 
invariant higher-dimensional subspaces. Using these de- 
generate eigen-states as a quantum memory, quantum in- 
formation storage of photonic states can be implemented 
up to a non-abelian holonomy. 
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FIG. 1. Four-level atom interacting with a quantum probe 
field (with coupling constant g, frequency w, and the detuning 
A p ) and two classic control fields (with frequency Vk , coupling 
Rabi frequency Qk, and the detuning Ak = u> a k — Vk, k = 1, 2). 
When Sk (= Ak — A p ) are very tiny, the system satisfies the 
near 3-photon resonance EIT condition. 

Our system consists of N identical 4-level atoms 
[13,14], where all the atoms are coupled to two single- 
mode classical control fields and a quantum probe field 
as shown in Fig. 1. The atomic levels are labelled 
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as the ground state \b), the excited state \a), and the 
meta-stable states \k) (k — 1,2). The atomic transition 
| a) — > |6), with energy level difference w Q b=w a — w&, is 
coupled to the probe field of frequency u> (= w a h — A p ) 
with the coupling coefficient g; and the atomic transi- 
tion | a) — > |fc) (fc = 1,2), with energy level difference 
w a fc, is driven by the classical control field of frequency 
Uk (= ^ak — A/j) with Rabi- frequency Ctk(t). 

In the present work we consider the case of 6k (= 
Afc — A p ) being very small, that is, those three fields have 
almost the same detuning with respect to the upper level 
| a). In view of the physical intuition, each meta-stable 
state with its relevant control field would constitute a 
near two-photon resonance EIT if another meta-stable 
state and its relevant control field do not exist. With the 
case of two-photon resonance EIT [17-19] (where the con- 
trol and probe fields have the same detuning) in mind, we 
would refer to our case of A p ~ A^ as a near "3-photon 
resonance" EIT. 

Under the rotating wave approximation the interaction 
Hamiltonian can be written as (let h = 1) [8] 

Hi = A p S + gy/NaA f + fi x exp[i0i (t)]Tj 1} 

+r2 2 exp[^ 2 (t)]T} 2) +h.c, (1) 

where 
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are symmetrized collective atomic operators. Here = 
\n)jj{v\ denotes the flip operator of the j-th atom from 
state \v)j to \fi)j (fi,v = a, 6, 1,2); and a the cre- 
ation and annihilation operators of quantum probe field 
respectively; and (pk(t) — Skt. The coupling coefficients g 
and Oi,2 are real and assumed to be identical for different 
atoms in the ensemble. A similar effective Hamiltonian 
was given in Ref . [8] for the case of an " atomic crystal" , in 
terms of quasi-spin-wave type collective atomic operators 
and a hidden dynamical symmetry was there discovered. 
The symmetrized operators (2) are just a special instance 
of the quasi-spin-wave operators discussed in [8] . 

Let us first consider a similar dynamical symmetry in 
the low excitation regime of the atomic ensemble where 
most of N atoms stay in the ground state \b) and N — > oo. 

(k) (k) 

It is obvious that T_ and Tj. (k = 1, 2) generate two 
mutually commuting SU(2) subalgebras of SU(3) [20]. 
To form a closed algebra containing St/ (3) and {A, A^}, 
we need to introduce two additional collective operators 

1 N 

^ = "/f fc =!' 2 ( 3 ) 



along with their hermitian conjugates. These opera- 
tors have the non- vanishing commutation relations Ck — 
[A,T^ k) ] and [C k ,T { _ k) ] = A (k = 1,2). As a special case 
of quasi-spin wave excitation with zero varying phases, 
the above three mode symmetrized excitations defined 
by A and C1.2 behave as three independent bosons. In- 
deed one can check that the operators (3), in the large N 
limit with low excitation, satisfy the bosonic commuta- 
tion relations [20]. The commutation relations between 
the St/ (3) algebra and the Heiscnbcrg-Weyl algebra h 
generated by A, A\ Ck, and CI imply that the dynamical 
symmetry of evolution governed by Hi can be described 
by the semi-direct product algebra SU (3)®/i. 

Based on this hidden dynamical symmetry of the in- 
teraction Hamiltonian, we can introduce a dark-state po- 
lariton operator 

D = a cos 6-C sin 6, (4) 

where 

C = C\ exp[i<f)i(t)] cos k + C2 exp[i<p2(t)] sinn (5) 

is a coherent mixing of two collective atomic excitations 
C\ and C2, and k = arctan^, 9 = arctan 2 ^^, = 
^/fif-i-fif. In terms of a new operator by 

T+ = exp[^i (t)] cos k + r{ 2) cxp[^ 2 (/)] sin k, (6) 

we can then rewrite the interaction Hamiltonian as 

Hj = A p S + gVNaA* + ttT+ + /i.e. . (7) 

Since [C, T_] = A and [A, T + ) = C, one can readily verify 
that [D, Hj] = 0. To generate the full eigen-space of H] 
with zero eigenvalue, i.e., the dark-polariton manifold, we 
need to consider another dark-state polariton operator 
complementary to D: 

E = C2 exp[z</5> 2 (£)] cos k — C\ exp[i<pi(t)] sinK. (8) 

It is worthwhile to point out that E satisfies the bosonic 
commutation relation as well and it is independent of D 
since [E, £)t] = 0. Moreover we have [E, Hj] = by con- 
struction. Our instantaneous quantum storage subspace 
V(t) is given by the linear span of the following family of 
instantaneous dark-states, i.e., the eigen-states of Hj(t) 
with vanishing eigenvalues 

\D m , n (t)) = -^=D^ n \0), (9) 

where |0) = |0) p (g) |b) = |0) p ® \b, 6, b) represents the 
ground state of the total coupled system with each atom 
being in the ground state \b) and the quantum probe 
field being in the vacuum state |0) p . It is easy to prove 
that any other dark-state polariton operator can be ex- 
pressed as a linear combination of D and E. Notice 
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that one can introduce the bright-state polariton oper- 
ator: B = as'm9 + C cos 9, which can be used to gener- 
ate eigen-states involving the excited state \a). However, 
as shown in Ref. [8], under the adiabatic manipulations, 
these states will not get coupled to the above constructed 
dark states (9). Of course, the states obtained by apply- 
ing D are not absolutely dark since the excited state can 
spontaneously decay. 

To study the geometric quantum information storage 
based on the above zero-eigenvalue dark states (9), we 
now consider the quantum evolution in V(t) caused by the 
adiabatic change of the parameters Clh (t) . The adiabatic 
condition is here given by [21,22] 
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for k = 1,2. Let us consider a state vector 
!*(*)) = E m .nC mn (t)\D m . n (t)) belonging to V(t). A 
straightforward calculation gives the matrix equation 
[21,22] for the coefficients c mn (t): d t C(t) = K(t)C(t), 
where the vector C(t) of coefficients and the connec- 
tion matrix K(t) are respectively defined by C(t) = 
(coo(*),coi(i), ...;ci (t),cn(t), ...) T , and K(t) m>m > tn>n > = 
-(D m >, n > {t)\d t D m . n {t)) (to, to', n, n' = 0, 1, 2, . . .). 

The quantum storage space V(t) is, in the considered 
limit, an infinite dimensional one. Thus in general it is 
difficult to write down the relevant connection matrix 
K(t) explicitly. On the other hand the adiabatic quan- 
tum evolution in V(t) can be reduced, i.e., this space 
splits into dynamically invariant finite-dimensional sec- 
tors. Let us explain this point now. 

We first observe that the following dynamical com- 
mutation relations hold fooif) '■= [D,dtD*] = 
-isin 2 6»((5iCos 2 K + S 2 sm 2 K ); ,f ED {t) := [E,d t D^] = 
— ksm9 + i(#2 — Si) sin # cos k sin k; foE{t) = 

l-fED(t)}* := [D,dtEfl]; f EE (t) := [E,d t E^] = 
— i(<5isin 2 K + (^2 cos 2 k). Using these relations for I' > 
to > and I > n > 0, we obtain 

(A'-m,m(t)|3 t A-n,n(i)) 

= Si> t id m>n [(l - m)f DD {t) + mf E E(t)} (11) 
+ 6i',i$ m ,n-i + l)(l - m)f DE (t) 

t 01' ,10 m,n+l 

yjm(l to + T)f E D (*) ■ 

Now it is clear from this expression that the total space 
V(t) can be decomposed into a direct sum of sub-spaces: 
V(t) - ®T=oMt), where Vj(t) = span{|D,-m,rr»(t)>|m = 
0,1,...,/} has dimension (I + 1). Notice that each 
Vi(t) is an invariant sub-space under the adiabatic ma- 
nipulation, i.e., if |$*(0)) e Vi(t), then \$i(t)) = 

J2 m c m( t )\ D l-m,m('t)) € V t (t). The restricted dynam- 
ics in Vi (t) is governed by the reduced dynamic equation 
dtCi(t) = K;(t)C;(t), where the sub-coefficient vector 
Ci(t) and the reduced connection matrix K;(£) are re- 
spectively given by d(t) = (%\t), cf \t)) T , 



and Ki(t) = (-(A-m,m(t)|9tA-n,n(<)))m,n=0,l,2,...,J- 

The solution Q(t) = W;(i) Q(0) formally determines 
the non-abelian holonomy W;(t) = Texp[/ Kj(t)dt], 
where T is the time-ordering operator. This non-abelian 
holonomy is non-diagonal and thus can mix different in- 
stantaneous eigen-states |D;_ TOim (t)) (m = 0, ...,l) in- 
ducing in this way a truly non-abelian gauge structure. 

In the following discussion, we consider the simplified 
model related to the above system as shown in Fig. 1: 
<5i,2 = 0. Such a system has only two controllable pa- 
rameters fii^ and can be readily realized experimentally. 
Mathematically the sub-connection can be simplified as 

K,(t) = HA- m ,m(t)|&A-n,n(t)>)m,T.=0,l,2,...,J (12) 



k sin 9 K 

'(0) 



(0) 
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where Kj ' is a constant matrix whose (to, n) entry is 



S m ,n-W ijn + l)(l - to) - 5 m , n+ iy/m{l - to + 1). (13) 

In this case the time-ordering becomes irrelevant and the 
non-abelian holonomy can be explicitly computed. In 
fact we have 



W,(t) = exp[^(t)K} 



(0)i 



(14) 



where 4>{t) — J ksm9dt. By noticing that is pro- 

portional to (E^D — h.c.) restricted Vi(t), we find that 
this non-abelian holonomy can be rather easily cast in 
a diagonal form by introducing a new instantaneous ba- 
sis. Let us introduce the new set of dark-state polariton 
operators 

D' = ^=(iD + E),E' = ^=(-iD + E), (15) 
and the associated dark states \D[_ n n {t)) = D J ^ |0). 

' ^/ ((— n).n. 

A straightforward calculation then gives a diagonal con- 
nection matrix K;(t) = — iks'm9dia,g (I, I — 2, —I) and 
the corresponding holonomy 



W{(t)=diag(e- 



-il<p(t) -i{l-2)4>(t) 



(16) 



Generally in EIT based quantum information storage 
protocols, the Rabi frequencies Qi ; 2 of the two classical 
control fields are initially set to a very large value com- 
pared to gy/N and then decreased independently and adi- 
abatically (e.g., as shown in Fig. 2). Thus 9{t = 0) — > 
and D(0) -> a. The initial state |$(0)) = X); c^WIOp® 
|b) can be written as |*(0)) = E*, m {0)\Dl_ m>m (0)), 
relative to the new basis \D[_ m m {t)) with the coefficients 



(iV2) l y/Tn\(l-m)\' 



(17) 



and [<&(*)) = Ej.m^WID^W), where c'«\t) = 
cxp[— i(l — 2m)(j)(t)]cm\o). When f2i j2 become negli- 
gible compared to gVN at time r, 9{t) — ► 7r/2 and 
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D(t) — > — Ci cos k — C2 sin k. This means that the quan- 
tum information, initially encoded in photonic states, has 
been transferred and written to atomic collective excita- 
tions. This accomplishes the quantum information stor- 
age protocol. 
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FIG. 2. Cyclic evolution of the parameters f2i,2. At time 
t = r, fii, 2 < g\/N; at time t = or t = T, fh,2 > gVN . 

In order to recover the stored information one needs 
to drive adiabatically the system parameters fii^ along 
a cyclic evolution such that at the time T the condition 
^1,2 gVN is satisfied in order to guarantee 9(T) — > 
(see the Fig. 2). At the intermediate times t £ (0, T) 
quantum information is encoded in a combination of pho- 
tonic and atomic collective excitations. In general, if one 
wants to recover exactly the initial state after that the 
adiabatic loop has been completed, she has to perform a 
unitary transformation to get rid of the effect of the non- 
abelian Berry phase factor. In particular, for a cyclic 
evolution of the parameters fii^ if 

<f>(T) = [ ksm6dt = 2j'K (18) 
Jo 

(J is an integer), it then follows that dm (T) = Cm (0). 
In this case the system state at the final time T coincides 
with the initial state |$(0)). 

We are now in the position, before concluding, to make 
a few comments on the relations between the results pre- 
sented in this paper and the general holonomic approach 
to quantum information processing [12,14]. In that ap- 
proach information is encoded in degenerate eigenstates 
of a parametric family of Hamiltonians, and in the generic 
case universal quantum computation [23] can be achieved 
by resorting to non-abelian holonomies only [12]. By re- 
garding the non-trivial holonomy one gets after an adia- 
batic loop as a designed quantum state transformation, 
rather than something one wants to get rid of, it should 
be then evident that the EIT-based scheme here dis- 
cussed represents an instance of such general strategy. 
For example the one exciton space Vi can encode one 
qubit: |0) := E^\0), |1) := D't|0). In this language the 
transformation (16) is nothing but a diagonal phase-shift 
[23]. In order to get non-diagonal single-qubit opera- 
tions one would have to relax the condition <5i = <5 2 = 0. 
Encoding many-qubit states and enacting a controllable 
geometric coupling between them - as required for realiz- 
ing universal computations - along with the robustness of 
the scheme against the various sources of errors is a more 
complex problem that calls for further investigations. 
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